A new Skyrme interaction with improved spin-isospin properties 
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A correct determination of the spin-isospin properties of the nuclear effective interaction should 
lead, among other improvements, to an accurate description of the Gamow- Teller Resonance (GTR). 
These nuclear excitations impact on a variety of physical processes: from the response in charge- 
exchange reactions of nuclei naturally present in the Earth, to the description of the stellar nu- 
cleosynthesis, and of the pre-supernova explosion core-collapse evolution of massive stars in the 
Universe. A reliable description of the GTR provides also stringent tests for neutrinoless double- 
/3 decay calculations. We present a new Skyrme interaction as accurate as previous forces in the 
description of finite nuclei and of uniform matter properties around saturation density, and that 



account well for the GTR in 48 Ca, u 
Dipole Resonance in 90 Zr and 208 Pb. 

PACS numbers: 21.60.Jz, 24.30. Cz 



Zr and Pb, the Isobaric Analog Resonance and the Spin 
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The Skyrme Hartrcc-Fock (HF) approach is one of 
the successful techniques for the study of the ground 
state properties of nuclei and, if supplemented by a 
proper description of nuclear superfluidity (e.g., within 
the Hartree-Fock-Bogoliubov scheme), it can be applied 
throughout the whole periodic table [1[ . The small ampli- 
tude limit of time-dependent HF calculations, or Random 
Phase Approximation (RPA), has allowed to describe 
many kinds of nuclear collective motion Q ■ The versatil- 
ity of the Skyrme ansatz allows its use in more elaborated 
theoretical frameworks that include higher-order nuclear 
correlations - like the Generator Coordinate Method Q , 
or the Particle- Vibration Coupling approach 

Despite the existence of drawbacks and open issues, 
the Skyrme-HF approach enables an effective description 
of the nuclear many-body problem in terms of a local 
energy density functional. Problems concerning specific 
terms of this functional need to be understood and even- 
tually solved (cf. also Ref. [f|). One of these problems, 
and the focus of the present work, is the accurate deter- 
mination of the spin-isospin properties of the Skyrme ef- 
fective interaction and of the associated functional. Such 
a determination should lead to accurate predictions of the 
properties of GTR, that are among the clearest manifes- 
tations of nuclear collective motion Q. Gamow- Teller 
(GT) transitions determine weak-interaction rates be- 
tween /p-shell nuclei that play an essential role in the 
core-collapse dynamics of massive stars leading to super- 
nova explosion 0, S| (in this neutron-rich environment, 
neutrino-induced nucleosynthesis may take place via GT 
processes Q). Accurate GT matrix elements are neces- 
sary for the study of double-/? decay [l(|, and may be 
useful in the calibration of detectors ai ming to measure 
electron- neutrinos coming from the Sun 

The earliest attempt to give a quantitative description 
of the GTR data was provided by the Skyrme SGII in- 



teraction 12|. However, two component spin-orbit con- 
tributions to the nuclear Hamiltonian density - Eq. (6.1) 
in Ref. (l3j - were not introduced. Later on, using the 
full Skyrme Hamiltonian density, a functional suitable 
for the predictions of finite nuclei and charge-exchange 
resonances was proposed: namely SkO' [l^. Relativistic 
mean-field and relativistic HF calculations of the GTR 

, 1||. In a rc- 
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have also become available meanwhile 
cent and detailed study on the GTR and the spin-isospin 
Landau-Migdal parameter (G' ) using several Skyrme sets 
[13], Bender et al. have concluded that this spin-isospin 
coupling is not the only important quantity in determin- 
ing the strength and excitation energy of the GTR in 
nuclei. Actually, the authors state that spin-orbit split- 
tings together with the residual spin-isospin interaction 
influences the above mentioned quantities. 

The GT transition strength (i?QT±) is mediated by 
the operator X^iLi cr (*) r ±(*)j where A is the mass num- 
ber and er and r are the spin and isospin Pauli matrices, 
respectively. The contributions to i?QT± come from nu- 
cleon transitions that change the spin and isospin of the 
parent quantum state and the residual interaction be- 
tween them is repulsive. The dominant transitions will 
be those between spin-orbit partner levels. In this re- 
spect, most of the Skyrme interactions overestimate the 
experimental spin-orbit splittings in heavy nuclei (l8j . 

Experimentally, the GTR exhausts only about 60- 
70% of the well known Ikeda Sum Rule (ISR) given by 
J[R GT -{E) - R GT +(E)]dE = 3{N - Z). To explain 
this well-known quenching problem, it has been proposed 
that the effects of the second-order configuration mixing, 
namely 2-particle 2-hole (2p-2h) correlations, or of the 
coupling with the A— hole excitation, have to be taken 
into account. The experimental analysis of 90 Zr [l9| 
seems to indicate that most of the quenching (around 
2/3) has to be attributed to 2p-2h coupling while the 
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TABLE I: Data and pseudo-data, Oi, adopted errors for the 
fit AOi, as well as partial and total number of data points 
and contributions to the \ 2 ■ 



Oi AO, 


v' 2 

Apartial 


^data 


Ref. 


B 1.00 MeV 


32.45 


5 


[20] 


r c 0.01 fm 


13.38 


4 


[21J 


AE so 0.04xOi 


19.02 


2 


[22] 


e„(p) 0.20x0; 


12.60 


11 


[23] 


x 2 


77.45 


/ 22 


= 3.52 



role played by the A isobar is much smaller. 

In our work, we present a new non-relativistic func- 
tional of the Skyrme type that include the central tensor 
terms (J 2 -terms) and two spin-orbit parameters. It is 
named SAMi for S'-kyrme ^4-izu Mi-lano. The new func- 
tional is as accurate as previous Skyrme models in the 
description of uniform nuclear matter properties around 
saturation and of masses and charge radii of double- 
magic nuclei. It is also precise in the description of the 
Giant Monopole Resonance (GMR), the Giant Dipolc 
Resonance (GDR) in 208 Pb and the GTR, the Isobaric 
Analog Resonance (IAR) and the Spin Dipolc Resonance 
(SDR) in medium and heavy mass nuclei. 

To this end, we have carefully chosen the set of fitted 
data and pseudo-data insp ired by the protocol used to 
build SLy interactions [13]: (i) the binding energies (B) 
of 40 ' 48 Ca, 90 Zr, 132 Sn and 208 Pb and the charge radii 
(r c ) of 40 < 48 Ca, 90 Zr and 208 Pb which allow us to de- 
termine the saturation energy (eoo), density (poo) and, 
to a good approximation the incompressibility (K^) of 
symmetric nuclear matter; (ii) the spin-orbit splittings 
(AE so ) of the lg and 2/ proton levels in 90 Zr and 208 Pb, 
respectively, that are well determined due to the flexibil- 
ity of our two-component spin-orbit potential; and (iii) 
the Landau-Migdal parameters Go and G' - associated 
with the spin and spin-isospin particle-hole (p-h) inter- 
action (see their definition in Ref. [3]) that are fixed at 
the values 0.15 and 0.35, respectively, at saturation den- 
sity. This features allow the new SAMi interaction to give 
an adequate description of spin-isospin resonances. In 
the literature, an empirical determination of the Landau- 
Migdal parameters can be found in Ref. t 25] but we do 
not use such values as pseudo-data in our fit. The reason 
is that the extraction of these values is based on single- 
particle energies obtained with a Woods-Saxon potential. 
In our case, we use HF energies associated with a differ- 
ent effective mass. However, we took inspiration form the 
empirical indications that suggest G' > Go > 0. This 
is not a very common feature within available Skyrme 
forces (see Fig. 1 in Ref. [24[). Therefore, we imposed 
that Go is larger than and G' is larger than 0.25, and 
we have tried to explore the optimum values that do not 
spoil the global fit; finally, (iv) pseudo-data correspond- 
ing to more fundamental microscopic calculations of the 
energy per particle of uniform neutron matter (e„) at 



baryon density p between 0.07 fm -3 and 0.4 fm -3 that 
have been helpful in driving the magnitude (J) and slope 
(L) of the nuclear symmetry energy at normal densities 
towards reasonable values. Table U provides references 
for these data and pseudo-data with the corresponding 
adopted errors, partial contributions to the x 2 , and the 
number of data points (ridata) used in the fit. The main 
differences between the SLy [l3| and the present pro- 
tocol are the fitting of the above mentioned spin-orbit 
splittings, the fact that we fix the spin and spin-isospin 
Landau-Migdal parameters, and the larger adopted er- 
rors for the equation of state of pure neutron matter. 
This protocol is justified by the fact that pseudo-data are 
used as a guide and, therefore, it should not impact on 
the fitted interaction more than experimental data. The 
minimization of the \ 2 has been performed by means of 
a variable metric method included in the MINUIT pack- 
age of Ref. [H. 

The parameters and saturation properties of the new 
interaction are shown in Table [TTJ The estimation of the 
standard deviation j2?| associated to each of them is also 
displayed. In what follows, the new SAMi functional is 
compared to available experimental data and other the- 
oretical predictions for ground and excited state prop- 
erties. First of all, we show in Fig. Q] the results for 
the symmetric and pure neutron matter Equations of 
State (EoS) as predicted by the benchmark microscopic 
calculations used in the fit [23], three state-of-the-art 
Brueckner-Hartree-Fock (BHF) calculations 28|-3(|, the 
SAMi functional, and SLy5 [13j - also fitted to repro- 
duce the neutron matter EoS of Ref. [23|. The agree- 
ment of the SAMi functional with these calculations of 
nuclear matter based on realistic nucleon-nucleon (NN) 
forces is remarkable. The deviation of the SAMi EoS of 
pure neutron matter from the fitted microscopic curve 
(red circles) is essentially due to the relatively large error 
(20%) adopted in the x 2 definition in order not to spoil 
the additional constraints set on the isovector channel 
of the effective interaction. In addition, we have checked 



TABLE II: SAMi parameter set and saturation properties (see 
text) with the estimated standard deviations [271 ] inside paren- 
thesis (referred to the last digits). 





value(cr) 




value (cr ) 


to 


-1877.75(75) MeV fm a 


Poc 


0.159(1) fm -3 


tl 


475.6(1.4) MeVfm 5 




-15.93(9) MeV 


t-> 


-85.2(1.0) MeV fm 5 


mis 


0.6752(3) 


t 3 


10219.6(7.6) MeV fm 3+3a 


rrirv 


0.664(13) 




0.320(16) 


J 


28(1) MeV 


J'l 


-0.532(70) 


L 


44(7) MeV 


£C 2 


-0.014(15) 


K x 


245(1) MeV 


x 3 


0.688(30) 


Go 


0.15 (fixed) 


Wo 


137(11) 


Go 


0.35 (fixed) 


Wo 


42(22) 






a 


0.25614(37) 
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FIG. 1: Neutron and symmetric matter EoS as predicted by 
the HF SAMi (dashed line) and SLy5 (solid line) interactions 
and by the benchmark microscopic calculations of Ref. [23|] 
(circles). State-of-the-art BHF calculations are shown by di- 
amonds [2S|, triangles [2!| and squares [30l |. 

that the SAMi EoS is stable against spin and spin-isospin 
instabilities [24| up to a baryon density of 4.1 p^ and 
5.3poo, respectively, i.e., well above the region important 
for the description of finite nuclei and enough for the 
study of uniform neutron-rich matter in neutron stars. 
Furthermore, we are aware that particle-number projec- 
tion techniques lead to instabilities when functionals with 
non-integer power of the density are employed (3lj . At 
the same time, the adopted density dependence (p a with 
a smaller than one) seems to be the only way to have 
reasonable values of the nuclear incomprcssibility and of 
the GMR energies within the Skyrmc functional. As a fu- 
ture perspective, all practitioners of local, Skyrme type, 
EDFs may need to deal with the problem of reproducing 
reasonable monopole energies on the one side and mak- 
ing particle number restoration doable on the other side. 
This is beyond the purpose of the current work. 

We display in Table Mil the SAMi results for binding 
energies, charge radii and proton spin-orbit splittings of 
all measured doubly-magic spherical nuclei. The descrip- 



TABLE III: Experimental data [2jj422] and SAMi results for 
the binding energies (B), charge radii (r c ), and proton spin- 
orbit splittings [AiJso (level)] on doubly- magic spherical nu- 
clei. 



El. 


A 


B 




r c 


' c 


AE so AE^f 


(level) 






[MeV] 


[MeV] 


[fm] 


[fm] 


[MeV] 


[MeV] 




O 


16 


130.48 


127.62 


2.75 


2.70 


5.26 


6.18 


(IP) 


Ca 


40 


347.08 


342.05 


3.47 


3.48 


5.78 


6.75 


(Id) 




48 


415.61 


415.99 


3.51 


3.47 


4.75 


5.30 


(Id) 


Ni 


56 


469.73 


483.99 


3.80 




4.69 


6.82 


(1/) 


Ni 


68 


593.19 


590.41 


3.92 










Zr 


90 


781.26 


783.89 


4.27 4.27 


6.45 


5.56 


(Iff) 


Sn 


100 


811.66 


824.79 


4.50 




4.60 


7.00 


(Iff) 




132 


1103.09 


1102.85 


4.73 




5.13 


6.51 


(lh) 


Pb 


208 


1636.61 


1636.43 


5.50 


5.50 


1.88 


2.02 


(2/) 



tions of B and r c are accurate within a 1% and 0.6%, re- 
spectively, and the proton spin-orbit splittings of the dif- 
ferent single particle states with high angular momenta 
are accurate within 15%. Such an accuracy will be deci- 
sive for an accurate characterization of the GTR energies 
171 and clearly improves the results obtained with SGII 

3- 

In Fig. [2] we test the performance of the SAMi in- 
teraction for the description of the strength distribution 
(calculated within RPA .321) in the cases of the GMR 
and GDR in 208 Pb. The results are compared with ex- 
perimental data and with the predictions of SLy5 [i~3j ]. 
The operators used in the GMR and GDR cases are, 

respectively, Y^=i r i and Z / A En=i r « - N / A Ep=i V 
The experimental centroid energy of the GMR has al- 
lowed to constrain the nuclear matter incompressibility 
at the value = 240 ± 20 MeV, by means of an anal- 
ysis of a large set of Skyrme interactions [33|. Within 



the same spirit, the experimental data on the GDR has 
allowed to determine the nuclear symmetry energy at a 
sub-saturation density S(p = 0.1 fm" 3 ) = 24.1±0.8 MeV 
341 ] . The SAMi interaction predicts compatible values, 
namely K x = 245 MeV and S(p = 0.1 fur 3 ) = 22 
MeV. Consistently, the giant resonance centroid energy 
predicted by SAMi agrees well with the experimental 
findings: £ c s AMi(G MR ) = 14.48 MeV should be com- 
pared with £° xp (GMR) = 14.24 ±0.11 MeV [ex- 
hausting both almost 100% of the Energy Weighted 
Sum Rule (EWSR) between E x = 8 - 22 MeV], and 
£f AMi (GDR) = 13.95 MeV should be compared with 
£^ xp (GDR) = 13.25 ± 0.10 MeV (exhausting both 
around 95% of the EWSR between E x = 9 - 20 MeV). 

The strength distributions of the GTR are displayed 
in Fig. |21 HF+RPA results obtained with the forces 
SAMi, SGII [13, SLy5 [H and SkO' 0] are compared 
with experiment. In the upper panel of Fig. [3l we show 
the experimental data of Ref. [37[ as well as the predic- 
tion of the SAMi, SGII and SkO' functionals for 48 Ca. 
In this case the SLy5 result is not shown because RPA 
produce instabilities. The nice agreement of the SAMi 
prediction in the excitation energy, E° xp = 10.5 MeV 
and £SAMi = W 2 MeV for the high -energy peak and 
Rf p = 3.0 MeV and £^ AMi = 2.0 MeV for the low- 
energy peak, and the % of the ISR exhausted by the 
main peak between 5 and 17 MeV, around 46% in the 
experiment and 71% in the calculation, is noticeable (in 
keeping with the fact that RPA does not include 2p-2h 
couplings). The prediction of the SAMi interaction in 
the case of 90 Zr (middle panel of Fig. [3]) is even bet- 
ter than in the case of 48 Ca. Despite of the accuracy 
of SGII, SLy5 and SkO' in describing other properties 
of nuclei, they do not perform as well as our new pro- 
posed functional. The excitation energy and % of the 
ISR exhausted by the high- and low-energy peaks in the 
experimental data (lil . |38| (in the calculation done with 
the SAMi functional) are, respectively, E° xp = 15.8 ±0.5 



4 



12 



> 

S 9 



I 

2 3 



i — 1 — i — 1 — ri — 1 — i — 1 — r 

— SLy5 i Exp. I 

- - SAMi 




10 12 14 16 18 
E [MeV] 



1 1 1 1 1 1 
; — SLy5 


I 1 . 
Exp. '. 


" - - SAMi 


t 208 

Pb : 


- /' 

- I 
I 


I 
11 1 
\i 


J 
I 

J 


\ 

\\ , '. 




! \\ \ 


Jl 


! V-A v - 
i \y \ \ 


// 


! \\ 




1 . • 1 . • . 


8 12 


16 



- 20 



- 15 



25 



> 



10 



Q 
O 



FIG. 2: Strength function at the rel- 
evant excitation energies in 208 Pb as 
predicted by SLy5 QJ] and the SAMi 
interaction for GMR (left panel) and 
GDR (right panel). A Lorentzian 
smearing parameter equal to 1 MeV 
is used. Experimental data for 
the centroid energies is also shown: 
£ C (GMR) = 14.24±0.11 MeV [31 and 
£ C (GDR) = 13.25 ± 0.10 MeV [3j|. 
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FIG. 3: GT strength distributions in 48 Ca (upper panel), 90 Zr 
(middle panel) and 208 Pb (lower panel) as measured in the 
experiment 0, and predicted by SLy5, SkO', SGII 

and SAMi forces. 

MeV and 57% (E$ Am = 15.5 MeV and 70%) between 
12 and 30 MeV and E^ = 9.0 ± 0.5 MeV and 12% 
(£SAMi = 7 8 MeV and 27%) between 3 and 12 MeV. 
In the lower panel and with unprecedented accuracy in 
HF+RPA calculations, the SAMi functional perfectly re- 
produces the excitation energy of the experimental GTR 
in 208 Pb fH]: E c **> = 19.2 ± 0.2 MeV and E% AMi = 19.3 
MeV. We also compare our results with the predictions 
of SGII, SLy5 and SkO' that fail in the description of 
the GTR in 208 Pb. It is important to notice that, op- 



posite to SLy5, the spin-orbit parameters (Wo and Wq) 
arc not fixed to be equal in the SAMi and SkO' inter- 
actions. Note also that, Go and G' were fixed to be 
0.011 and 0.503 in SGII interaction 01 together with 
A"oo=215 MeV and J=26.8 MeV which give reasonable 
descriptions for other resonances but predict the GT ex- 
citation energies in 208 Pb at slightly higher values (21.2 
MeV) - partly due to a larger G' . 

In order to assess the robustness of the SAMi func- 
tional in the description of other charge-exchange re- 
actions, we have analyzed the Isobaric Analog Reso- 
nance (IAR) in 90 Zr and 208 Pb. In this nuclear exci- 
tation, modclized by the operator X^=i r ±W; the nu- 
cleon transitions change the isospin of the parent quan- 
tum state. We have found that the experimental value 
for the IAR in both, 90 Zr and 208 Pb is very well repro- 
duced by the SAMi functional - within a 1.5% discrep- 
ancy with respect to the experiment [1H - and the sum 
rule, J[Ria-{E) - R IA +(E)]dE = (N - Z), is perfectly 
exhausted in our calculations. 

Furthermore, experimental and theoretical studies on 
the Spin Dipole Resonance (SDR) have been recently re- 
vitalized due to its connection, via a sum rule, to the neu- 
tron skin thickness of nuclei (Ar„ p ) 4l| and, therefore, to 
the density dependence of the nuclear symmetry energy 
0,S3- For these reasons, we also present the SAMi pre- 
dictions for this important charge- exchange excitation in 
90 Zr (Fig. [H) and 208 Pb (Fig. [5]) as compared with the ex- 
periment [40, EH . The operator used for the RPA calcu- 
lations is Yyf=i Em T ±(*) r f i Y L(fi) ® cr(i)] JM and, as it 
is shown in both figures, it connects single particle states 
differing by a total angular momentum: J n = 0~,1 ~~ 
and 2". The sum rule, J[R SD -(E) - R SD +(E)]dE = 
^r(iV(r 2 } — Z(r 2 }), is completely exhausted in our cal- 
culations, 99.99% in the case of 90 Zr and 100% in the 
case of 208 Pb. The experimental (calculated) value for 
the sum rule is 148 ± 12 fm 2 [IH (150 fm 2 ) for the case 
of 90 Zr. The total and multipole decomposition of the 
experimental (calculated) value for the integral of 
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FIG. 4: Spin Dipole strength distributions in 90 Zr as a func- 
tion of the excitation energy E x in the r_ channel (upper 
panel) and r+ channel (lower panel) measured in the experi- 
ment [4lJ and predicted by SAMi. Multipole decomposition 
is also shown. A Lorentzian smearing parameter equal to 2 
MeV is used. 

i? SD - in the case of 208 Pb are 1004± 24 f m 2 ( 12 24 fm 2 ), 
107±? fm 2 (158 fm 2 ) for J 71 = 0", 450±^ fm 2 (423 fm 2 ) 
for J 71 = 1" and 447±^ fm 2 (643 fm 2 ) for J 71 = 2". 
Finally, it is noticeable the overall agreement between 
experiment and SAMi predictions when the strength dis- 
tributions as a function of the excitation energies shown 
in Figs. @] and [5] arc compared. 

The neutron skin thickness of medium and heavy nuclei 
is known to be strongly correlated with the isospin prop- 
erties of the nuclear effective interaction [42|, [43| ■ A re- 
cent study [44| shows that the Ar np in 208 Pb derived from 
different hadronic probes agrees in a value of 0.18±0.03 
fm. The SAMi interaction predicts 0.15 fm, compatible 
within the estimated error bars. In addition, a recent the- 
oretical study [HI has allowed to predict 0.17 ± 0.02 fm. 
Recently, the PREx collaboration has reported a value of 
0.33^0^8 f m f° r the- same observable measured via par- 



ity violating elastic electron scattering [46j, |47J. If this 
value is confirmed with high accuracy, a deep revision of 
current nuclear models will be necessary. In the case of 
90 Zr, an analysis of the charge exchange Spin-Dipole res- 
onance has allowed to extract Ar„ p ( 90 Zr) = 0.07 ± 0.04 
fm j4lj . in perfect accordance with our predicted value 
of 0.07 fm: an additional proof of the improvement in 
the description of the spin and isospin channels of the 
nuclear effective interaction provided by SAMi. Finally, 
the neutron skin in 48 Ca is predicted by our model to be 
0.17 fm very close to the theoretical value 0.176 ±0.018 
fm reported in 45 1. 

In summary, we have successfully determined a new 
Skyrme energy density functional which accounts for 
the most relevant quantities in order to improve the 
description of charge-exchange nuclear resonances, i.e., 
the hierarchy and positive values of the spin and spin- 
isospin Landau-Migdal parameters and the proton spin- 
orbit splittings of different high angular momenta single- 
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FIG. 5: SDR strength distributions for 208 Pb in the r_ 
channel from experiment [ZTJ] and SAMi calculations. To- 
tal and multipole decomposition of the SDR strength are 
shown: total (upper panel), ,V — CP (middle-upper panel), 
J n = 1~ (middle-lower panel) and J n = 2~ (lower panel). A 
Lorentzian smearing parameter equal to 2 MeV is used. 

particle levels. As a proof, the GTR in 48 Ca and the 
GTR, IAR and SDR in 90 Zr and 208 Pb arc predicted 
with high accuracy by SAMi without deteriorating the 
description of other nuclear observables and, therefore, 
promising its wide applicability in nuclear physics and 
astrophysics. 
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